Proper orthogonal decomposition of solar photospheric motions 
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Abstract 

The spatio-temporal dynamics of the solar photosphere is studied by performing a Proper Or- 
thogonal Decomposition (POD) of line of sight velocity fields computed from high resolution data 
coming from the MDI/SOHO instrument. Using this technique, we are able to identify and char- 
acterize the different dynamical regimes acting in the system. Low frequency oscillations, with 
frequencies in the range 20-130 //Hz, dominate the most energetic POD modes (excluding solar 
rotation), and are characterized by spatial patterns with typical scales of about 3 Mm. Patterns 
with larger typical scales of ~ 10 Mm, are associated to p-modes oscillations at frequencies of about 
3000 /iHz. 



PACS numbers: 96.60.Mz; 92.60.Ek; 96.60.Ly; 89.75.-k 



The solar photosphere is an interesting example of a system exhibiting complex spatio- 
temporal behavior, a quite common feature in a wide range of systems far from equilibrium. 
The early research in pattern formation focused on the presence of simple periodic struc- 
tures, while the main questions currently addressed concern regimes characterized by higher 
complexity, that is, patterns that are more irregular in space and time. This is often related 
to the occurrence of intermediate states between order and turbulence j])]. High resolution 
images of the solar photosphere show that the whole surface of the Sun is covered by a 
cellular pattern, the so called solar granulation, consisting in bright cells (granules) , with 
roughly 10^ km size, surrounded by a dark background of intergranular lanes cf. P, 0]. 
This structure represents the manifestation of the complex convective dynamics, at high 
Rayleigh numbers, occurring in the subphotospheric layers jsl |^. Other patterns are also 
observed in the photosphere, namely the supergranulation 3 x 10^ km size) P, 0] and 
giant cells of size about 10^ Km, perhaps the result of global convection over the whole 
solar convective depth In the solar photosphere the convective dynamics coexist with 
the contribution of solar oscillations. The interior of the Sun, because of the gravitational 
field and density stratification, behaves like a resonant cavity and supports the excitation of 
global oscillations . Depending on the restoring force, two main types of waves can be 

excited, namely acoustic high frequencies p-modes [sj, 1^ (in the range 1000-5000 /iHz jl^ ) 
and gravitational low frequency g-modes P, ^ (in the range 1-200 /xHz) . The Fourier k-u 
spectra of line of sight photosperic velocity fields show ridges corresponding to the discrete 
p-modes ^ (see figE))- On the other hand, discrete frequencies in the low-frequency range 
(g-modes) are not recognized, but rather the power observed in this range is spread over a 
continuum, commonly attributed to the solar turbulent convection. Indications for the pres- 
ence of 160 min discrete oscillations within a continuum have been reported |8j, followed by 



contradictory evidences 



12j . Despite these attempts, no unambiguous mode identification 



and classification has been established so far. 

In the present paper, we use the Proper Orthogonal Decomposition (POD) as a power- 
ful tool to investigate the dynamics of stochastic spatio-temporal fields. In astrophysical 
contexts, POD has been recently used to analyze the spatio-temporal dynamics of the solar 
cycle Introduced in the context of turbulence [l^, the POD decomposes a field u{r,t) 
as u{r,t) = Yl'jLo^ji^)'^ ji^) (-'-)' eigenfunctions being constructed by maximizing the 
average projection of the field onto constrained to the unitary norm. Averaging leads to 
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an optimization problem that (^(r, t), ^(r', t))\l/(r')dr' = A\l/(r) (2), where 

Q represents the spatial domain and brackets represent time averages. The integral equa- 
tion (2) provides the eigenfunctions and a countable, infinite set of ordered eigenvalues 

> -^j+i) each representing the kinetic energy of the j-th mode. Thus, POD builds up the 
basis functions, which are not given a priori, but rather obtained from observations. The 
time coefficients aj{t) are then computed from the projection of the data on the correspond- 
ing basis functions ^E'j(r) so that the sum (1), when truncated to terms, contains the 
largest possible energy with respect to any other linear decomposition of the same trunca- 
tion order. The POD basis functions are optimal in turbulence with respect to the classical 
Fourier analysis, where the basis functions are not proper eigenfunctions of the signal. This 
method is particularly appropriate when analyzing complex physical systems, where differ- 
ent dynamical regimes coexist. POD allows to identify these regimes and to characterize 
their energetics and their spatial structure. We show that POD is able to capture the main 
aspects of the spatio-temporal dynamics of the solar photosphere. 

The line of sight velocity fields u{x, y, t) {x, y being the coordinates on the surface of the 
Sun) used in this work have been obtained from images acquired by the Michelson Doppler 
Imager (MDI) instrument mounted on the SOHO spacecraft The image size is 695 x 695 
pixels, with a spatial resolution of about 0.6 arcsec (1 arcsec ~ 725 km on the surface). The 
time series spans a time interval T = 886 min, images being sampled every minute. The 
top panel of Figure Q shows one of the snapshots of the dataset, together with its spatial 
spectrum. These data have been already studied in the past with Fourier techniques 
and the k-u power spectrum of the velocity signal is shown in figure El 

The POD of the velocity field u{x,y,t) yields a set of eigenfunctions y) and coeffi- 
cients aj(t), as well as the sequence of eigenvalues A-,- (j = 0, 1, . . . , 886, sorted in decreasing 
energetic content). Most of the energy is associated with the first POD mode, accounting 
for the line of sight component of solar rotation. As this is neither the expression of solar 
turbulence nor oscillations, this first POD mode will be ignored throughout this paper. The 
rest of the energy is decreasingly shared by the following 885 modes, so that, for example, 
only 4% is associated to the second (j = 1) POD mode. In laboratory turbulent fiows, POD 
attribute almost 90% of total energy to the typical large-scale coherent structures, confined 



to j < 2 



17| . In our case, 90% of energy (excluding the rotation) is contained in the first 



140 modes. This indicates the absence of dominating large-scale structures and the presence 
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of some turbulent dynamics related to nonlinear interactions among different modes and 
structures at all scales. 

We turn now the attention to the eigenfunctions, namely the spatial patterns, found by 
POD, and reported in figure Q for three modes. The power spectra |\E'j(A;)^| obtained from 
Fourier transform of the eigenfunctions are also reported in figure ^ for the same 

three modes, as a function of the wavevector k = (k^ + kyY^'^. It can be observed that the 
patterns are rather disordered, with a broad spectrum. Moreover, modes can be separated 
in three different groups, accordingly to grain size, as is shown in figure H The 1 < j < 11, 
most energetic modes display a pattern with very fine structures, recalling the photospheric 
granulation, with a broad spectrum indicating the complex nature of the POD modes, and 
show a center-limb modulation. Conversely, eigenfunctions of the modes j = 12 and j = 13 
present a coarser pattern, and their (still broad) spectra show a number of ridges. The 
further lower energy modes cannot be precisely classified, and seem to present a mixture 
of the previous characteristics, both in grain size and in spectral shape. The spatial pat- 
tern associated to each mode, although complicated, can be quantitatively characterized by 
computing the integral scale length Lj = iV'j(^) 
the energy containing scale of classical turbulence 
scale for the first ten fine grained modes L ~ 3 Mm, while for the coarse grained modes 
L ^ 10 Mm. 



k ^dk/ dk which represents 
3]. This allows to estimate the typical 



FIG. 1: Left column, from top to bottom: one snapshot of the original data, and the POD 
eigenfunctions ^j{x,y) in the plane {x,y) for the three modes j = 8, j = 12 and j = 49. Each 
image covers about 300 x 300 Mm^ on the solar surface. Right column: the wave vector spectra 
|^'j(A:)p of the corresponding images. On the second and third spectra, the wave vector spectra 
obtained from the k-u; cut, respectively at u; = 76 /xHz and uj = 3254 /iHz, are superimposed 
(dotted lines, in arbitrary units). 



FIG. 2: The k-cu spectrum of the dataset. 



The eigenfunctions described above are associated with the corresponding coefficients, 
accounting for the temporal evolution of each mode. In figure El we report, for the same 
modes represented in figure ^ the time behaviour of the coefficients aj{t), together with 
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their Fourier spectra |aj(z/)p. As can be seen, for all the eigenfunctions the spectra of the 
time coefficients present two peaks, located in two well defined and separated frequency 
ranges. The fine grained eigenfunctions (e. g. j = 8) are associated with time coefficients 
dominated by low-frequency oscillations. For the coarse grained patterns (e. g. j = 12) 
the high frequency peak becomes of the same order, and even slightly higher than the low 
frequency one. For the intermediate cases (e. g. j = 49), the amplitudes of the two peaks 
are again of the same order, but with weak prevalence of the low frequencies. For each 
POD mode, the high frequency peaks can be identified from the spectra, lying in the range 
3250-3550 /iHz. In a similar way, for the low frequency oscillations, we find frequencies in 
the range 20-127 /iHz. 

FIG. 3: Left column: time evolution of POD coefficients aj{t) (m/s) for the same three modes as 
in figure n Right column: the corresponding frequency spectra |aj(z^)p {rn? /s^). 

We try now to understand the results observed with POD. Let us first focus on the high- 
frequency modes. Both the informations on the involved spatial scales and the measured 
frequencies indicate that solar p-mode contributions are identified by POD. Indeed, the fre- 
quencies we obtain (for instance, v = 3405 uHz for ?' = 12, z/ = 3367 uHz for ?' = 13) are 

nnn 

compatible with the p- modes observed using Fourier techniques p|, lUJ, , and have been 
predicted by helioseismological models cf. j9|. The spatial pattern of about 10 Mm, which 
POD associates with high frequency modes, is in agreement with the horizontal coherence 
length attributed to solar p-modes. In order to confirm this, we compared the k-u classical 
results with the POD spectra. Choosing, for example, the POD mode j = 12, we performed 
a cut of the k-uj spectrum at the measured peak frequency (3254 yuHz). We thus obtain a 
wavelength spectrum, reported on figure ^ that can be compared with the j = 12 eigen- 
function spectrum. As can be seen, the spikes observed in the POD spectrum qualitatively 
correspond to the typical ridges structure of the k-u spectrum. 

Concerning the low-frequency dominated modes, it is interesting to note that the mea- 
sured frequency range, mentioned above, is compatible with the theoretical results on solar 
gravitational modes 0, 0]. However, because of the resolution limit, related to the finite 
length of the time series, the low frequency oscillations detected by POD can not be un- 
ambiguosly attributed to discrete modes, rather than to a continuum. The evidence of 
center-limb modulation shows that such modes are mainly associated with the contribution 



of horizontal velocities. To confirm this, we tried to apply the POD technique to a reduced 
square of about 100 x 100 pixels taken at the center of each image, where horizontal com- 
ponents of the velocity are weaker. In this case, the low-frequency oscillations are no longer 
observed (see figure 0]). As for the p-modes, the k-uj cut can be performed at the peak fre- 
quency of, for example, mode j = 8 (76 /iHz). The comparison with the POD eigenfunction 
wavelength spectrum is shown in figure ^ and again reveals a qualitative agreement between 
the two techniques. The characteristic integral scale computed for the k-uj cut is also in 
agreement with the corresponding POD integral scale. Since the POD basis functions are 
not chosen a priori, the association between small-scale eigenfunctions and low-frequencies 
is reliable. In order to test this, we reduced the resolution of the data by performing a 
spatial 50 x 50 pixels running average of the fields before applying POD. In this case, the 
low-frequency peaks in the coefficient spectra are completely lost (see figure E)). 

FIG. 4: Upper panels: the POD eigenfunction ^q{x, y) (left, the center-limb modulation is evident) 
and the Fourier spectra (right) obtained from two reduced square boxes of size 100 x 100 pixels, 
taken at the center (dotted) and at the border (full) of the image (see left panel). Bottom panels: 
50 X 50 pixels running averaged data at a given time (left), and the j = 9 POD mode wave vector 
spectrum (right) for both the smoothed field (full) and the original field (dotted). 

In conclusion, we presented the first application of POD on high resolution solar photo- 
spheric velocity fields, which represents an example of convective turbulence in high Rayleigh 
numbers natural fiuids. POD is able to capture the main energetic and spatial features of the 
solar photosphere. In particular, the dynamical processes of interest are automatically sep- 
arated from the most energetic, but uninteresting mode witch captures solar rotation. Two 
main oscillatory processes, weU separated in frequency, are detected in ah the other modes: 
high frequency oscillations in the range 3250-3550 /iHz and low frequency oscillations in 
the range 20-130 /iHz. The high frequency waves are the well known acoustic p-modes and 
their properties, as obtained from POD, are in agreement with previous results based on 
Fourier techniques. On the other hand, low frequency oscillations, with frequencies com- 
patible with those expected from the theory of solar g-modes, prevail in the most energetic 
POD modes (excluding solar rotation), which are characterized by spatial eigenfunctions 
with typical scales of ~ 3 Mm. The clear association between low frequency oscillations and 
small spatial scales, which are close to the solar granulation, is the most interesting and most 
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surprising result provided by our analysis. This suggests the presence of a strong coupling 
between low frequency modes and the turbulent convection, which could have implications 

n 

on the inversion of time-distance measurements of sound waves |21], also used to infer the 
properties of the turbulent convection below the surface. Such coupling could be related to 
the dishomogeneities and small scale structures arising from the highly nonlinear dynamics 
of the convective layers. This point needs to be investigated in future theoretical studies. 
A further improvement of this analysis can be expected from its application to longer time 
series. 
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